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ABSTRACT 

Context. Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and 
powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree 
of the wind mass-loss depends on the initial metallicity of WR stars. 

Aims. Following our comprehensive studies of WR stars in the Milky Way, M 31, and the LMC, we derive stellar parameters and mass- 
loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity 
environment on the mass loss and evolution of WR stars. 

Methods. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The 
physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations. 

Results. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar 
temperatures exceeding 75 kK, while their luminosities range from 10^ ^ to 10®' Lq. The WN stars in the SMC exhibit on average 
lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M 31, and the LMC. 

Conclusions. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear 
dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by 
radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high 
temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface 
hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past. 

Key words. Stars: Wolf-Rayet - Magellanic Clouds - Stars: early type - Stars: atmospheres - Stars: winds, outflows - Stars: 
mass-loss 


1. Introduction 

The Small Magellanic Cloud (SMC) is a unique laboratory for 
studying massive stars in a low-metallicity environment. With 
a distance of 60.2 kpc, corresponding to a distance modulus of 
18.9 mag (|Westerlund||1997 [Harries et al.||2003| |Hilditch et al.| 


2005 


Keller & Wood|2006 , it is close enough to allow for de¬ 


tailed spectroscopy of its massive star population. Moreover, the 
interstellar extinction in the direction and within the SMC is 
on average quite low (e.g.. Grieve & Madore||1986 Schwer- 


ing & Israel||1991 Bessell||1991| Schlegel et al. 1998| [Larsen 


et aH2000[ [Haschke et al.|2011[ l. The metallicity of SMC stars 


is in general subsolar, but the individual values are strongly age 
dependent ( [Piatti|20TT] [2012| l. Even in comparison to the metal 
deficient LMC, the metal content of SMC stars is still on average 
a factor of 2.5 lower at all stellar ages ( Piatti & Geisler[2013[ ). 

Stellar winds of hot stars are assumed to be driven by the 
transfer of momentum from photons to the matter of the stellar 
atmosphere by scattering and absorption in spectral lines (for a 
review see, e.g., [Puls et al.|2008] l. Eor this reason, the strength 
of the winds and the associated mass-loss is a function of the 
opacity and should thus depend on the metal content of the stellar 
atmosphere (e.g., Kudritzki et al. [1W7[ Puls et al!]|2000[ Vink 


et al. 2000 2001[ l. Eor Wolf-Rayet (WR) stars, the dependence of 


by [Grafener & Hamann[ ( [2008[ l and Vink & de Koter[ ( [2005| l. 
Since the SMC has a subsolar metallicity, considerably weaker 
radiation-driven stellar winds are expected in the SMC than in 
the Galaxy. 

It is well known that in the SMC, the WR stars of the nitro¬ 
gen sequence (WN stars) exhibit emission line spectra with rel¬ 
atively weak and narrow emission lines (e.g., [Conti et 3T 1989| l. 
[Crowther & Hadfield| ( [2006| l found that the line luminosities of 
the WN stars in the SMC are a factor of four to five lower than 
their LMC counterparts. Moreover, distinct absorption lines can 
be recognized in eight of the twelve WR stars known in the SMC. 
These characteristics were either interpreted as a sign of binarity 
( [Massey et al.|2003] l or attributed to photospheric lines that are 
inherent to the SMC WN stars because of their relatively weak 


stellar winds (Conti et al. 


1989 


[et al.[ ( [2003| l and Eoellmi (2004 1 showed that these absorption 


Massey et al. 20031. Eoellmi 


lines are blue-shifted but that they exhibit a constant radial ve¬ 
locity, concluding that these lines must form in the wind of these 
objects. Later on, reduced mass-loss rates of the SMC WN stars 
(compared to their LMC and MW counterparts) were supported 
by the findings of Crowther ( 2006[l, Nugis et al. (2007 1 , and Bo- 


[nanos et al.| ( [2010 1 . 


the mass-loss rate on the metallicity was theoretically predicted 


The WR population of the SMC is small but exceptional. 
It comprises only twelve objects. In contrast more than 140 WR 
stars are known in the LMC, and more than 630 in the Galaxy. In 
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the SMC, almost all WR stars belong to the WN sequence. The 
only exception is the binary system SMC AB 8 whose primary 
belongs to the rare class of WR stars showing prominent oxygen 
lines in their spectra (WO type). Wolf-Rayet stars of the carbon 
sequence (WC stars) are completely absent in the SMC. This is 
in strong contrast to the Galactic WR population where the num¬ 
ber of WC and WN stars is approximately equal. For compari¬ 
son, the WN/WC ratio in the LMC is roughly four. These dif¬ 
ferences in the specific WR populations are generally attributed 
to the different metallicities of these galaxies. The same trend of 
an increasing WN/WC ratio with decreasing metallicity is also 
found in other Local Group Galaxies (|L6pez-Sanchez & Esteban 

|20T0t . 

In addition to the WN/WC ratio also the distribution of the 
subtypes within each sequence varies between galaxies of dif¬ 
ferent metallicities. Using stellar atmosphere models, |Crowther| 
& HadfieTd](|2006|l showed that a down scaling of the mass-loss 


rate at lower metallicity has a direct influence on the spectral 
appearance of WR stars, resulting in weaker emission line spec¬ 
tra and earlier WN subtypes for significantly subsolar metallic¬ 
ity objects. Indeed, almost all SMC WN stars exhibit an early 
spectral subtype (WNE: WN2-WN5), the only exceptions being 
SMC AB 4 and the secondary in the SMC AB 5 binary system. In 
contrast, the ratio of early- to late-type WN stars (WNL: WN6- 
WNll) in the MW is almost unity, while the majority of the 
LMC WN stars belong to the WNE category. 

Despite these characteristic differences between WR popula¬ 
tions of individual galaxies, the impact of the metallicity on the 
strength of WR winds is not yet well established by means of 
spectral analyses. |Crowther| ( |2006| l and |Nugis et al.| ( |2007| l em¬ 
pirically studied the wind mass-loss as a function of the metal¬ 
licity based on a subset of the WR stars known in the Galaxy, 
LMC, and SMC. [Vink & de Koterj ( |20051 l derived a theoretical 
prediction for the dependence of the mass-loss rate of WNL stars 
on the metallicity by means of Monte-Carlo simulations. 

Empirical and theoretical studies indicate a power-law de¬ 
pendence of the mass-loss rate M on the metallicity Z with an 
exponent in the range of 0.8-1.0. However, a comprehensive 
evaluation of the metallicity effect on WR winds, considering 
a large number of objects from the Local Group Galaxies is still 
missing but imperative. The present study addresses this impor¬ 
tant problem, building on our recent analyses of WN stars in the 
Galaxy (|Hamann et al.|2006|l, M31 (|Sander et al.|2014|l, and the 
LMC dHaimch et al.|2014^ 

The paper is structured as follows: In Sect.j^we briefly char¬ 
acterize the sample investigated in this work before presenting 
the data our analysis is based on in Sect.|^ We then describe 
the Potsdam Wolf-Rayet (PoWR) atmosphere models and the 
procedure of the analysis in Sect.|^ The results are compiled in 
Sect.|^ In Sect. 1^ we discuss the evolutionary status of the WN 
stars in the SMC. Section]?] gives a summary and our general 
conclusions. The appendices encompass additional tables (Ap¬ 
pendix]^, comments on the individual objects (Appendix]B|i, 
and spectral fits of all stars in our sample (Appendix]C|l. 


2. The sample 

The number of known WR stars in the SMC has increased since 


the first comprehensive compilation by Azzopardi & Breysacher 
from eight to twelve ([Morgan et al.||199T]~ Massey & 


Duffy] 2001] Massey et al. 20b3| l. We identify each star by 


( |1979) l and which is one of the suggested names by the SIM- 
BAD database. Alias names of the objects are included in the 
comments on individual objects (Appendix]B]l. 

In the present study, we concentrate on the putatively sin¬ 
gle WN stars in the SMC. |Poellmi et al.| ( |2003| l and Poellmi| 
( |20041 l carefully studied the binary status of the SMC WN stars 
by means of radial-velocity shifts, line profile variations, pho¬ 
tometric measurements, and X-ray data. According to these au¬ 
thors, SMC AB 1,2,4, 9, 10, 11, and 12 are probably single stars. 
We kept SMC AB 9 in our sample although ]Poellmi et al.| ( |2003) 
found weak indications that this object might actually be a bi¬ 
nary. Thus, the sample investigated in this work comprises seven 
objects. The spatial positions of these stars in the SMC are shown 
in Pig.]^ A detailed analysis of the known binary WR systems 
in the SMC will be the subject of a forthcoming paper (Shenar 
et al., in prep). 


3. Observational data 

Eor most of our targets, high-resolution spectra taken with the 
Ultraviolet and Visual Echelle Spectrograph (UVES) mounted 
on Unit Telescope number 2 (UT2) of the Very Large Telescope 
(VLT) are available. The fully reduced UVES data in a wave¬ 
length range of 3960^960 A and 5040-6010 A were retrieved 
from the ESO archive for advanced data products. Eor details on 
the data reduction we refer the reader to the ESO phase 3 data 
release websit^ We rectified the spectra by means of the spec¬ 
tral continua which were identified by a careful adjustment of a 
low-order polynomial to the observed data. 

Eor those objects without available UVES sp ectra, we used 
low-re solution spectra (6.7-2.8 A) obtained by 


Eoellmi et al. 


a 

running number, following the name scheme SMC AB # which 
is based on the work conducted by [Azzopardi & Breysacher] 


( 2003 j l in the spectral range of 3700-6830 A. Although these 
datasets were primarily intended for radial velocity studies, they 
are well suited for our analysis owing to their high signal-to- 
noise ratio (S/N). Detailed information on the instrumentation 
used and the data reduction can be found in jEoellmi et al.| ( ]2003| l. 
These spectra were normalized and co-added by the former au¬ 
thors. Eor three objects in our sample, we also used optical spec¬ 
tra observed by jTorres-Dodgen & Masse3^ ( |1988| l. In contrast to 
the spectra described above, these datasets are flux-calibrated. 
While the S/N of these low-resolution spectra is relatively mod¬ 
est, they are only used to construct the observed spectral energy 
distribution (SED). 

Additionally, archival spectra taken with the Ear Ultravio¬ 
let Spectroscopic Explorer (EUSE) and the International Ultra¬ 
violet Explorer (lUE) were retrieved from the MAST archive. 
High-resolution EUSE spectroscopy in the EUV (960-1190 A) 
is available for four objects in our sample. To improve the S/N 
ratio, these spectra were rebinned to 0.1 A, resulting in a spectral 
resolution which is still sufficient to resolve the interstellar ab¬ 
sorption lines. Low-dispersion lUE spectra (resolution of ~ 6 A) 
in the range of 1200-3300 A were downloaded for three stars. 
Since the lUE and EUSE data are flux-calibrated, these spectra 
are rectified by the reddened model continuum. A compilation 
of the spectra used in the spectral analysis of each sample star 
can be found in Table lA.2l 

In addition to flux calibrated spectra, we used narrow-band 
magnitudes given by Crowther & Hadfieldj ( |2006| l and [Torres-j 
[Dodgen & Massed ( |1988| l to construct the observed SED in the 
optical range. Photometric measurements in the I band are taken 


* https://WWW.eso.org/observing/dfo/quality/PHOENIX/ 
UVES_ECH/processing.html 
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R. Hainich et al.: Wolf-Rayet stars in the Small Magellanic Cloud 



Fig. 1. The position of the putatively single WN stars in the SMC. The stars are identified by their number, following the scheme SMCAB# 
(see text for details). The background image shows the Out nebular emission in the SMC. The image was retrieved from the Magellanic Cloud 
Emission-Line Survey (MCELS, [Smith et al.|2005^ . 


from [Monet et al.| ( |2003| l and [Bonanos et al.] ( |2010| l. In the near 
infrared we used photometry (J,H,Ks) from the Two Micron 
All Sky Survey (2MASS, Cutri et ar||2012| l. Furthermore, we 
applied Spitzer-IRAC magnitudes published by [Bonanos et al.| 
(|2010|l and the SAGE EMC and SMC IRAC Source Catalog 
(IPAC 2009). Eor SMC AB 2, 4, & 11 we also incorporated pho¬ 
tometry from the WISE all sky survey ( [Cutri & et al.|2012] l. Ad¬ 
ditional photometry was used for SMC AB 11 as described in 
Appendix[Bl 


4. The analysis 

Eor the spectral analysis we used synthetic spectra calculated 
with the PoWR model atmosphere code. By fitting the model 
spectra to the observational data, we derive the fundamental pa¬ 
rameters of the stars and their winds. 


4. 1. The models 

PoWR is a state-of-the-art code for expanding stellar atmo¬ 
spheres. Its basic assumption is a spherically symmetric, sta¬ 
tionary mass-loss. The models account for deviations from local 
thermodynamic equilibrium (non-LTE), wind inhomogeneities, 
and iron line blanketing. The radiative transfer equation is solved 
in the co-moving frame consistently with the rate equations for 
statistical equilibrium, while the energy conservation is also en¬ 


sured. Some details on the numerical methods and the assump¬ 
tions applied in the code can be found in |Grafener et al.| ( |2002) l 
and |Hamann & Grafener| ( |2003|[2004| l. 

The inner boundary of the model atmospheres is located at 
a Rosseland continuum optical depth of Tross = 20, which per 
definition corresponds to the stellar radius R*. The stellar tem¬ 
perature r* is related to this radius by the Stefan-Boltzmann law 


L = Ano-RlTt 


( 1 ) 


stellar temperature thus represents an effective temperature 
= 20 

OSS 

In the subsonic part of the stellar atmosphere, the velocity 
held is calculated self-consistently to fulfill the hydrodynamic 
equation ( [Sander et al.|2015] l. The velocity held in the supersonic 
part of the stellar wind is prescribed by a so-called /(-law ( [Castor 
[& Earners 1 1979 [[Pauldrach et al.|1986] l. Eor most of the models 
calculated in the context of this work, the exponent (3 is set to 


unity. Only in the case of SMCAB 10 is a double-yS law (Hillier 
[& Miller[1999 [Grafener & Hamann[2005[ ) in the form described 
by Todt et al~| ( 2015[ l used since it gives a slightly better ht (for 
details see Appendix[B|). 

Wind inhomogeneities are accounted for by means of the so- 
called mircoclumping approach, which assumes optically thin 
clumps ( |Hillier[1991[[Hamann & Koesterke[1998 1 that hll a vol¬ 
ume fraction /y, while the interclump medium is void. The den- 
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sity contrast between the clumps and an homogeneous model 
with the same mass-loss rate is given by the clumping factor 
D - /v~* ■ The models used here account for a depth dependency 
of D{r), where the clumping starts at the sonic point, increases 
outward, and reaches the full value of D at 30 % of the terminal 
velocity. 

The main wind parameters of a stellar atmosphere model 
are combined in the so-called transformed radius ( [Hamann] 
|& Koesterke|1998[ ) 


/ 

Voo 

I M'/D ' 

2500kms-i 

V 

/ 10-4Moyr-iJ 


( 2 ) 


relating the mass-loss rate M with the terminal wind velocity Voo, 
the stellar radius /?», and the clumping factor D. 

( |1989| l found that the equivalent width possess by 
in synthetic spectra is approximately independent from the spe¬ 
cific value of the luminosity and the individual values of the wind 
parameters which enters Rx as long as the transformed radius, the 
stellar temperature, and the chemical composition remain unal¬ 
tered. Therefore, models of WR stars can be easily scaled to dif¬ 
ferent luminosities. According to Eq. (|^, the mass-loss rate has 
to be scaled proportional to Under the additional constraint 
of a constant Voo also the line profiles are conserved. The scal¬ 
ing invariance for models of the same transformed radius can be 
understood from the fact that R^ is proportional to the ratio be¬ 
tween the emission measure of the wind and the stellar surface 
area. 

The model calculations rely on complex model atoms of H, 
He, N, C, and in some cases also O and P (see Table |A. l] for 
details). The iron group elements, comprising Fe, Sc, Ti, Cr, 
Mn, Co, and Ni, are treated in the “superlevel approach” (see 
[Grafener et al.|2002l l. 


Schmutz et al. 
emission hnes 


4.2. Abundances 


The atmospheres of WN stars exhibit burning products of the 
CNO cycle because of envelope stripping due to the stellar winds 
and internal mixing processes such as rotationally induced insta¬ 
bilities. Hence, the surface abundances of these objects are al¬ 
tered in comparison to that observed in unevolved stars such as 
B-type main sequence stars. Since the proton capture of is 
the slowest reaction in the CNO cycle, most of the initial carbon 
and oxygen is converted into nitrogen. According to |Schaerer| 
et al. (1993 i, the remaining carbon and oxygen is only 1/60 of 
the nitrogen mass fraction. 

By default, we thus assume a nitrogen abundance which is 
roughly the sum of the observed CNO abundances for B-type 
stars and Hit regions. For the initial composition of the SMC 
WN stars, we refer to stellar atmosphere analyses of B-type 
stars carried o ut by|Korn et al. ( |2000] l. Trundle et al. (2007 1 and 
Hunter et al.| ( 2007[ ). Additionally, we take the chemical com- 
position of H ii regions in the SMC derived by |Kurt & Dufour| 
( |1998| l into account. As described above, we calculate the chem¬ 
ical composition of our grid models from the averaged abun¬ 
dances obtained by these authors. Accordingly, we used a nitro¬ 
gen abundance of An = 1.5 ■ 10“^ (mass fraction) and a carbon 
abundance of Ac = 2.5 ■ 10“^. The iron abundances is taken 
from s tudies conducted by |Kom et al. 1 ( 12000) 1 and [Trundle et al.| 
(2007 1 , resulting in an averaged value of Ape = 3 ■ lO"'*. These 


basic model abundances, as well as the observed values obtained 
by the former authors, are listed in Table[T] 

For comparison, Table[T] also includes the abundances used 
by [Hamann et al.j ( |2006| l and [Hainich et al.j ( |2014| ) for the anal¬ 
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yses of the MW and LMC WN populations, respectively. The 
CNO abundance of the SMC grid models is almost a tenth of the 
solar value, while the iron abundance is approximately a fifth. 
Even compared to the LMC, the SMC CNO and iron abundances 
are lower by more than a factor of two. We note that the abun¬ 
dances for the SMC WN stars reflect the basic composition of 
the grid models. The nitrogen abundance for all analyzed stars 
was individually adjusted during the fitting procedure. For those 
stars which show oxygen lines in their spectra oxygen was also 
included in the model calculation. The individual abundances 
derived in our analyses can be found in Tables]^ and A. 3 We 
note that the values derived for the CNO and P abundances are 
uncertain by up to a factor of two. The error in the derived hy¬ 
drogen abundances is on the order of 5 %. 


4.3. Spectral fitting 

The spectral analysis is based on three grids of stellar atmo¬ 
sphere models. The free parameters of these grids are the stel¬ 
lar temperature T* and the transformed radius Rx, while the grids 
are distinguished by hydrogen abundances of 0.2, 0.4 and 0.6 
(mass fraction). For further details on the grid spacing we refer 
the reader to Hamann & Grafener| ( 2004] l. The model grids are 
publicly available on the PoWR homepage ^ Suitable model s are 
preselected by mean s of contour plot s ( Sander et al.|2012) l or a 
;^fy-fitting technique |Todt et ^ 2013 1 . Subsequently, individual 
models with refined stellar parameters and abundances are cal¬ 
culated for each object in our sample. 

A typical example for a normalized line fit is shown in Fig.|^ 
As mentioned in Sect.[T] distinct absorption lines are visible in 
the spectra of most SMC WN stars. These features are character¬ 
ized by asymmetric line profiles and a strong blue-shift ( jFoellmij 
et al.|2003 1 that does not vary in time, indicating the severe in¬ 


fluence of the stellar wind on these lines. With the exception of 
a few lines (He ii AA4200, 4542A in the spectra of SMC AB 1 
and 12), the absorption line component as well as the relatively 
weak emission-line component of the SMC WN-star spectra can 
be simultaneously reproduced by the models, illustrating that 
weak absorption lines do not necessarily imply binarity (see also 
[Hamann et al. [2006) 1. 

The stellar temperatures of the model atmospheres are ad¬ 
justed until the observed He i/He n and N iii/N iv/N v line ratios 
are adequately reproduced. For SMC AB 2 and 4, this procedure 
leaves a temperature uncertainty of approximately +3kK. For 
the other stars in our sample, the etTor in the temperature might 
be higher because the available spectra for these objects only 
possess lines of one ionization stage per element. The individual 
error for each object will be discussed in AppendixjB] 

The transformed radius is adjusted such that the observed 
strengths of the emission lines are reproduced by the model. The 
error margin of Rx is on the order of +0.05 dex, with the excep¬ 
tion of SMC AB 4, where it is twice as large. This translates to 
an uncertainty of the mass-loss rate of +0.075 dex (+0.15 dex for 
AB 4). According to Eqs. ([1]) and (|^, the uncertainty in the stel¬ 
lar temperature corresponds to an additional uncertainty in the 
mass-loss rate of +0.025 dex for the cooler stars in our sample, 
while it is twice as large for the hotter objects. The terminal ve¬ 
locity and the degree of clumping are derived from the width 
of the emission lines and the strength of the electron scattering 
wings ( [Hillier[[T991[ [Hamann & Koesterke[[1998) , respectively. 
The accuracies in the determination of the terminal velocity and 
the clumping factor are on the order of 10 % and 50 %, respec- 


www.astro.physik.uni-potsdam.de/PoWR.html 
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Table 1. Chemical composition (mass fractions in percent) 



Sun'* 

MW/M31 

WN'' 

EMC 

WN" 





B stars* 

SMC 

B stars" 

Hii^ 

WN« 

c 

0.237 

0.01 

0.0067 

0.02 

0.022 

0.014 

0.0025 

N 

0.069 

1.5 

0.40 

0.003 

0.018 

0.0037 

0.15 

O 

0.573 

- 

- 

0.13 

0.17 

0.12 

- 

SCNO 

0.88 

1.51 

0.41 

0.1323 

0.21 

0.1377 

0.1525 

Fe 

0.129 

0.14* 

0.07* 

0.035' 

0.027 

- 

0.03* 


Notes. 

lAsnlund et al. (l2009k 

([2007 

(e) 

Korn et al. ( 

2000pTl 

[Trund 

le et al.|(|20071( 


;t al. 1 2014^ as used in Hainich et al. i2014i Hunter et al. 
L this work including iron group elements mean value from 




X[A] 


Fig. 2. Normalized line spectrum of SMC AB 12 in comparison to the a PoWR model. The observation is shown by a thin (blue) solid line, whereas 
the best fitting model spectrum is overplotted by a thick (red) dashed line 


lively. These uncertainties contribute to the error margin of the 
mass-loss rate (M oc which is at most 0.2 dex. 


Simultaneously with the line fit, the simulated spectral en¬ 
ergy distribution (SED) is fitted to flux-calibrated spectra and 
photometry in order to determine the luminosity and the color 
excess (see Fig.|^ for an example). For this purpose the 
model flux is geometrically diluted according to a distance mod¬ 
ulus of 18.9 mag (|Westerlund|1997[[Harries et al.|2003[|Hilditch| 


et al.|2005] [Keller & Wood|2b06[ ). The reddening encompassed 

contributions from an internal SMC component and the Galactic 
foreground. For the latter we used the reddening law published 
by Seaton] ( 1979| l and a color excess of Et,^^ = 0.05 mag, which 
is the mean reddening derived by [Grieve & Mador^ ([1986[l, 


[Bessell[ ( 1991| l, [Schwering & Israel| ( [1991| l, and [Schlegel et al!] 
( [1998[ l. For the SMC component, the color excess is a free pa¬ 
rameter that is individually adjusted, while the reddening law is 
adopted from the work of[Gordon et al.[p003[). 


The error margin of the luminosity has a couple of contri¬ 
butions: the uncertainty in the color excess, which is relatively 
small for the SMC stars; the uncertainty in the temperature de¬ 
termination; the approximated treatment of the ISM extinction; 
and the error inherent to the photometry and the flux calibrated 
spectra. Since the flux of WN stars scales approximately linear 
with the stellar temperature in the optical and IR (Rayleigh-Jeans 
domain), the temperature uncertainty contributes 0.03 dex to the 
luminosity error of the WN 5-6 stars and about 0.1 dex to the er- 
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Fig. 3. The observed spectral energy distribution of SMCAB4, composed of flux calibrated spectra (FUSE, lUE, CTIO; blue thin lines) and 
photometry (blue boxes, labeled with magnitudes), compared to the model continuum (red dashed line) of the best fitting model. The model flux 
is corrected for interstellar extinction and accounts for the geometrical dilution due to the distance of the SMC. 


ror in the luminosity of the WN 3^ stars. The uncertainty in the 
SED fit results in a luminosity error of approximately 0.05 dex 
for those objects where both flux-calibrated spectra and pho¬ 
tometry is available, while it is twice as high for SED fits that 
are based on photometry alone. Therefore, the total error margin 
in the luminosity amounts to 0.1 dex for the cooler objects and 
0.2 dex for the hotter objects without flux calibrated spectra. 


5. Results 

5.1. Stellar parameters 

The basic stellar parameters derived for the SMC WN stars are 
compiled in Table|^ Eurthermore, the number of hydrogen and 
helium ionizing photons emitted by these objects and the associ¬ 


ated Zanstra temperatures are given in Table A.4 

In Eig.|^ we show the location of our program stars in the 
log T.-log/Jt-plane. Although the number of objects is rather 
small, a clear one-dimensional sequence is recognizable with 
SMC AB 4 being the only outlier. In relation to the other SMC 
stars, the position of the latter object in this diagram implies a 
significantly denser wind and higher mass-loss rate. Again with 
the exception of SMC AB 4, the surface hydrogen abundance of 
the SMC WN stars gradually decreases with increasing stellar 
temperature. Compared to their Galactic and EMC counterparts, 
the WN stars in the SMC have on average higher stellar temper¬ 
atures, which in turn reflects their subtype distribution. 

In Table]^ we also give the wind efficiency, wfiicfi is defined 
as tfie wind momentum divided by the momentum of the radia¬ 
tion field: 


'I = 


MvocC 


(3) 


Witfi their comparably weak winds, all WN stars in the SMC 
have wind efficiencies below unity. They thus clearly fall below 
the single scattering limit. In strong contrast, the early-subtype 
WN stars (WNE) have averaged wind efficiencies of 2.1 in tfie 
EMC ( [Hainicfi et al.|[20T4| ) and 5.9 in tfie WM ( |Hamann et al.| 
|2006| l. Hence they clearly exceed the so-called single scattering 
limit and enter the multi scattering domain, where photons need 
to be scattered rj times on average to drive the wind. 

Except for the weakness of their winds (see Sect. |53] l, the 
most intriguing property of all WN stars in the SMC is the oc¬ 
currence of hydrogen in their atmospheres. Eor SMCAB 1,2, 
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r,/kK 



log(r*/K) 


Fig. 4. The position of the analyzed SMC WN stars in the log T,-logRi- 
plane. The symbols reflect the WN subtype (see inlet). The hydrogen 
abundances is indicated in three steps by means of the filling color (pur¬ 
ple: Ah > 50 %, red: Ah = 30 - 50 %, orange: Ah < 30 %). 


and 4 this fact was already noticed by Conti et al. 


( 1989| l. Later, 


|Foellmi et al. ( |2003| l and |Foellmi| ( 2004| l showed that all puta¬ 
tively single WN stars in the SMC possess hydrogen at their 
surface. Our spectral analysis quantifies tfiese findings, yielding 
hydrogen abundances in the range of Ah = 0.2 - 0.55 (by mass). 
These values can be well understood in terms of weak stellar 
winds that were not able to completely remove the hydrogen rich 
layers of the stellar atmospheres of these stars. 

Early type WN stars that exhibit a high stellar tempera¬ 
ture, hydrogen at the stellar surface, and an inherent absorp¬ 
tion line component in their spectra can also be found in the 
EMC, although they are rare. Two examples are BAT99 25 and 
BAT99 74. Especially the latter exhibits striking similarities with 
SMC AB 1, which is, however, more luminous by a factor of two. 
For comparison, no early type WN stars with absorption fines 
and a significant amount of hydrogen are known in the MW, 
indicating that metallicity is decisive for the evolution of mas¬ 
sive stars (see Sect. |5.3.2| and [fill. In this context, BAT9925 and 
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BAT99 74 might represent objects with initial metallicities sig¬ 
nificantly below the standard LMC metallicity, explaining their 
spectral appearance and their unusually low mass-loss rates. 
This concept is corroborated by SMCAB4. In contrast to the 
other SMC WN stars, AB 4 exhibits spectral characteristics and 
a mass-loss rate that reflect rather typical LMC values, suggest¬ 
ing that the progenitor was more metal rich. 

Alternatively to the scenario stated above, interaction with 
a close stellar companion might explain the abundance pattern 
seen in some of these stars. Mass transfer via Roche-lobe over¬ 
flow is able to efficiently remove parts of of the hydrogen rich 
atmosphere of the donor star, resulting in a WR star with a reju¬ 
venated companion. This could explain, e.g., the relatively low 
hydrogen abundance derived for SMC AB 4 with respect to its 
late spectral type. However, the study conducted by |Foellmi et al.| 
( 2003| l revealed that the WN stars discussed in this paper do most 
likely not have a close companion, rendering it unlikely that 
these objects have been the donor stars in a binary system. Of 
course, more complex binary interactions might have happened. 
For example, some of the SMC WN stars could have been mass 
gainers in a binary interaction, while the system was disrupted by 
the supernova (SN) of a previous primary. In this case, however, 
it remains to be proven that the rejuvenated object can evolve 
into a WR star, which might be questionable at SMC metallici¬ 
ties. 


In Figs.j^ and we plot the absolute visual magnitudes 
(monochromatic colors as defined by Smith|1968) l of our sample 
stars versus their spectral subtypes (WN 2-WN 6). The hashed 
areas in these figures refer to the corresponding relations for the 
LMC (Hainich et al. 2014| l and MW (accounting only for the 
stars with known distances, |Hamann et al.| 2006) 1, respectively. 
Although the SMC stars appear to be slightly brighter in the vi¬ 
sual than the LMC stars, there is a fairly good agreement be¬ 
tween the results from the three galaxies, indicating a universal 
correlation between and spectral type. This is quite astonish¬ 
ing, given that the WNE stars from individual galaxies appear to 
be very different, e.g., most WN 2 to WN 4 stars in the LMC and 
MW are hydrogen-free, while all their SMC counterparts exhibit 
hydrogen at the stellar surface. 

For three stars in our sample (SMC AB 1, 2, 4), spectral anal¬ 
yses by means of stellar atmosphere models were conducted by 
Crowtlier](2000 1 and Martins et al. (|2009)l. In general, the basic 


stellar and wind parameters deduced in our analysis agree well 
with the results published by these authors. The only exception 
is SMCAB 1. For this object, we obtain a significantly higher 
stellar temperature and consequently also a higher luminosity 
(see Append ix[B]for deta ils) than the values derived by |Martins| 
et al. (|2009|l. Nugis et al. (20071 computed the luminosity for six 


stars in our sample by means of an adopted bolometric luminos¬ 
ity of 4.74 mag and bolometric corrections derived by |Nugis~&] 
Lamers (2000[l. In Table|A.5] we compare these estimates with 


the luminosity derived from our elaborated SED fits. Except for 
SMC AB 2, the luminosities obtained by |Nugis et al.| ( |2007| l are 
significantly lower. 


5.2. The Hertzsprung-Russell diagram 


The Hertzsprung-Russell diagram (HRD) of the single WR stars 
in the SMC is shown in Eig.[7| Eor comparison, the HRD posi- 


tions obtained by Hainich et al. P014 1 for the LMC WN stars 

and byjHamann et al. ( 

2006)1, jMartins et al. ( 

20081, Liermann 

et al.|(|20101l, and|Barnis 

ce et al.)()2008)l for the Galactic WN stars 


are also reproduced by gray and open symbols, respectively. Eor 



WN subtype 


Fig. 5. Absolute monochromatic u-magnitudes ( |Smith|[T968) of the 
SMC WN stars plotted over suhtype number. Mean magnitudes for each 
subtype are represented by thick blue symbols, while red dots refer to 
individual objects. The straight red line is the linear regression to the 
SMC stars. The hashed area represents the scatter in the LMC correla¬ 
tion published by|Hainich et al.|(2014^. 



Fig. 6. Same as Fig.[^ but the hashed area refers to the scatter in the 
Galactic relation. 


the SMC stars, the hydrogen abundance, as derived in this anal¬ 
ysis, is color-coded. 

With the exception of SMC AB 2, all SMC WNE stars are 
located on the hot side of the zero-age main sequence (ZAMS). 
These objects, while containing a considerable amount of hy¬ 
drogen in their stellar atmosphere, are located in a temperature 
range (> 79 kK), where almost none of their LMC and Galactic 
counterparts show a measurable amount of hydrogen. In addi¬ 
tion, these hot stars (except of SMCAB 10) exhibit luminosities 
of log (L/Lq) > 5.85 and thus occupy a parameter space where 
neither MW nor LMC WR stars can be found. Therefore, despite 
the low number statistics, one might conclude that the WNE stars 
in the SMC decisively differ from the Galactic and LMC popu¬ 
lation. In contrast, the HRD position and hydrogen abundance of 
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Table 2. Parameters of SMC WN stars 


SMCAB 

Subtype" 

T, 

[kK] 

logRt 

[Re] 

^co 

[km/s] 

Eb-v 

[mag] 

[mag] 

R. 

[Re] 

logM 

[Mo/yr] 

D 

logL 

[Tg] 

7 

M’’ 

[Me] 


Comments 

1 

WN3ha 

79 

1.37 

1700 

0.22 

-4.57 

5.7 

-5.58 

10 

6.07 

0.2 

75 

0.5 


2 

WN5ha 

47 

1.63 

900 

0.10 

-4.96 

9.1 

-5.75 

4 

5.57 

0.2 

43 

0.55 


4 

WN6h 

45 

1.30 

1000 

0.09 

-5.85 

13.0 

-5.18 

10 

5.78 

0.5 

38 

0.25 


9 

WN3ha 

100 

1.22 

1800 

0.09 

-3.94 

3.5 

-5.65 

10 

6.05 

0.2 

62 

0.35 


10 

WN3ha 

100 

1.18 

2000 

0.12 

-3.38 

2.2 

-5.64 

4 

5.65 

0.5 

37 

0.35 


11 

WN4h:a 

89 

1.35 

2200 

0.11 

-3.69 

3.5 

-5.56 

4 

5.85 

0.4 

50 

0.4 

composite 

12 

WN3ha 

112 

1.14 

1800 

0.14 

-3.34 

2.4 

-5.79 

10 

5.90 

0.2 

42 

0.2 



Notes. The subtype classifications follow Foellmi et al. 1 2003 i and Foellmi 1 2004 
mass-luminosity relation derived by|Grafener et al.H2011^ 


m 


Masses calculated from the luminosity, using the 


r»/kK 

200 150 too 60 40 20 10 



Fig. 7. The HRD of the single WN stars. The color-filled symbols repre¬ 
sent the WN stars in the SMC analyzed in the present paper. The filling 
color of these symbols reflects the derived atmospheric hydrogen abun¬ 
dance (see inlet). The gray filled symbols and the open symbols refer to 
the LMC WN stars and the Galactic WN stars, respectively. 


SMC AB 4, the only WNL star in the SMC, roughly corresponds 
to the WNL population in the MW and the LMC. 

None of the putatively single WN stars in the SMC has a 
luminosity of log (L/Lq) < 5.6, in sharp contrast to the WN pop¬ 
ulation in the MW and the LMC. This finding promotes a sig¬ 
nificantly higher minimum initial mass for the WR stars in the 
SMC, in line with theoretical expectations for its low metallic- 
ity (see Sect.|^. Interestingly, this effect is not seen between the 
MW and LMC ( [Hainich et al.|2014) l, although the difference in 
their metallicity content is also not negligible. 


MW, LMC, and SMC, showing that the weak emission lines of 
the SMC WN stars correspond to a reduced mass-loss rate. The 
present study confirms these findings. All our sample stars have 
mass-loss rates of about log(M/(MQ/yr)) a; -5.6, with the ex¬ 
ception of SMC AB 4, which displays a mass-loss rate roughly a 
factor of three higher than any other WN star in the SMC. 


In Table |A.5l we compare the mass-loss rates obtained from 
our detailed spectral line fits with the results presented by |Nugis| 
|et al.| ( |2007| l. These authors derived the mass-loss rates for six 
stars in our sample using the equivalent width of the He ii A 4686 
line and an empirically scaled formula. This comparison re¬ 
veals a reasonable agreement, exhibiting the largest deviations 
for SMC AB 1 and 2, where the values presented by |Nugis et al.| 
( 2007| l are larger by about 50 %. 


The mass-loss rates of our SMC WN sample and of the LMC 
WN stars ( [Hainich et al.||20T4| l are compared in Fig.|^ While 
the SMC stars are on average slightly more luminous than the 
bulk of the LMC stars, their mass-loss rates are, on average, 
lower by more than a factor of four. In Fig.w e also plot the 
mass-loss relations given by [Nugis & Lamersj ( |2000| l for WN 
stars (their Eq. 20) and for WR stars in general (their Eq. 22). 
These relations are based on mass-loss rates determined from 
radio emissions and emission line equivalent widths, accounting 
for the effect of clumping in both cases. While the universal WR 
calibration published by jNugis & Lamers] ( |2000| ) intrinsically ac¬ 
counts for a metallicity dependence, the WN relation does not. 
We therefore scaled their WN relation according to a metallicity 
dependence of M cx predicted by | Vink & de KotCT (2005 i. 

Both prescriptions are shown for surface hydrogen abundances 
of 40 % and 60 % (by mass) in Eig.|^ The individual values pre¬ 
dicted for our program stars can be found in Table|^ The third 
and fifth column of this Table give the factor between the mass- 
loss rates predicted and the results obtained from our detailed 
model fit. In general, the latter values are considerably lower 
than the prediction of [Nugis & Earner^ ( |2000| ), with SMC AB 2 
being the only exception. We note that the disagreement tends to 
be larger for lower surface hydrogen abundance. Such discrepan¬ 
cies might imply that the relation suggested by [Nugis & Lamersj 
(|2000|l is oversimplified. 


5.3.1. Wind-momentum luminosity relationship 


5.3. Stellar wind mass-loss 


As described above, the emission lines in SMC WR stars are 
systematically weaker than their Galactic and LMC counterparts 
( [Conti et al.|1989 Crowther & Hadfield[2006[ l. [Crowth^ ( [2006[ l 
and Nugis et al. (20071 studied a subset of the WN stars in the 


The “modified wind momentum” is a measure for the strength of 
the stellar wind ( [Kudritzki et al.[ 1995) [Puls et al.[1996[[Kudritzki| 
[et al.|199^ that is defined as 


T>mom = MVocRI^^ ■ 


(4) 


Article number, page 8 of 


30 







































































R. Hainich et al.: Wolf-Rayet stars in the Small Magellanic Cloud 



Fig. 8. The mass-loss rates as a function of the luminosity of the WN 
stars in the SMC as well as the LMC for the sake of comparison. The 
meaning of the symbols is the same as in Fig[7l The d ashed lines refer to 
the mass-loss relation published by |Nugis & Lamer^ ( |2000^ for Galactic 
WN s tars, which is scaled acc ording to a metallicity dependence of M oc 
Z® jVink & de Koter|2005 1. The straight lines are the general mass- 
loss prescription for WR stars given by |Nugis & Lamer^ ( |2000[ l. Both 
relations are plotted for two surface hydrogen abundances. 


Table 3. Mass-loss rates predicted by the relations from Nugis & 
|Lamers| ( |2000[ l 


SMC AB 

WN relation" 
logMN&L 
[Mo/yr] 

WR relation 
logMN&L 
[Mo/yr] 

1 

-4.93 

4.5 

-4.96 

4.2 

2 

-5.85 

0.8 

-5.68 

1.2 

4 

-5.01 

1.5 

-5.03 

1.4 

9 

-4.71 

8.7 

-4.79 

7.2 

10 

-5.36 

1.9 

-5.3 

2.2 

11 

-5.11 

2.8 

-5.11 

2.8 

12 

-4.56 

17.0 

-4.83 

9.1 


Notes. scale d according to a metallicity dependence of Z” “ i Wnk& 
de Koter|2Q05) ratio between the mass-loss rates predicted by Nugis 


& Lamers|(|2000^ and the values derived in our analysis 


The theory of radiation driven winds predicts a tight relation of 
the form L, the so-called wind-momentum luminosity 

relationship (WLR, Kudritzki et al.|1999] l. 

In Fig.|^ we plot the modified wind momenta versus the lu¬ 
minosities of the WN stars in the SMC and its LMC counter¬ 
parts ( [Hainich et al. 12014) . Furthermore, we show the empirical 
WLRs for OB-type stars derived by Mokiem et al.| ( p007] l, us¬ 
ing clumping corrected Dmom values for stars with Her emission. 
As expected from the theory of line driven winds, the WN stars 
in the SMC exhibit on average a lower wind momentum than 
their higher metallicity counterparts. This difference amounts 
to 0.45 dex for those objects with similar luminosities. Interest¬ 
ingly, the WN stars in the SMC posses modified wind momenta 



Fig. 9. The modified wind momentum Z)n,om as a function of the lu¬ 
minosity for the putatively single WN stars in the SMC (same sym¬ 
bols as in Fig.|^. For comparison the LMC WN stars (gray filled sym¬ 
bols) analyzed by jHainich et al.|j2014^ also shown. The blue and gray 
straight line is the clumping corrected wind luminosity relation for OB- 
type stars in the SMC and LMC, respectively, derived byjMokiem et al.j 

pboTi i. 


that are in the range of what is expected for equally luminous 
O-type stars. In contrast to the SMC stars, most of the LMC 
WN stars exhibit winds which are significantly stronger than the 
wind of equally luminous O-type stars (on average almost a fac¬ 
tor of six). Thus, the modified wind momentum supports the pic¬ 
ture that the wind of the WN stars in the SMC are considerably 
weaker than their counterparts in the LMC. 


5.3.2. Metallicity dependence 

As highlighted in Sect.[T] it is expected that the mass-loss rate of 
massive stars varies as a function of the metallicity Z. By apply¬ 
ing a modified CAK approach ( [Castor et al.|1975] Abbott|1982[ 
[Pauldrach et al. 19861 on O-type stars with Magellanic-Cloud 
metallicities, Kudritzki et al.[ ( [T987[ l derived a scaling relation of 
the form 


M~Z“ 


(5) 


with a - 0.5. Vink et al. ( 2000[ 200 Ij l obtained an exponent 
of a = 0.69 by means of a Monte-Carlo approach to the wind 
dynamics. This theoretical prediction is supported by empirical 
studies (e.g., Bouret et al.[2003][de Kote7| 2006 [Mokiein et al.[ 
2007|l. For WR-stars, the situation is less clear since only a lim¬ 


ited number of observational constraints ( [Crowther|2006[ |Nugis[ 
[et al.|2007| [Hainich et al.|2014[l and even fewer theoretical pre- 
dictions ( [Vink & de Koter|2005[ [Grafener & Hamann|2008[ ) are 

available. 

In recent years, a large number of WN stars were spectro- 


scopically analyzed in the MW (Hamann et al. 

2006 Martins 

[et al.|2008| Liermann et al. [20101 [Oskinova et al. 

20T35;inM31 


cilitating a comprehensive investigation of the relation between 
their wind mass-loss and metallicity. The averaged mass-loss 
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Fig. 10. The averaged mass-loss rates of the WN stars in the SMC 
(blue), the LMC (red), the MW (dark green), and the late type WN 
stars in M 31 (light green) are plotted versus the logarithmic metallcity 
of their host galaxies. The error bars refer to the standard deviation in 
each sample. The dashed line is a linear regression. 


Fig. 11. Averaged mass-loss rates of the WN stars in the SMC, LMC, 
MW, and the late type WN stars in M 31 plotted versus the iron abun¬ 
dance given in Table[T] The symbols are the same as in Fig.[T] The 
dashed line is a linear regression fit that is weighted by the number of 
objects analyzed in each galaxy. 


rates of the WN stars in the MW, the late type WN stars in M 31, 
the WN stars in the LMC, and the SMC WN stars analyzed in 
this work are plotted in Fig. 10 as a function of the metallicity of 
their host galaxy. For the MW, we assume that a metallicity of 
Zq = 0.014, derived for the solar neighborhood by |Przybilla et al.| 
( 2008| l, is valid for the whole Galaxy, neglecting the observed 
metallicity gradient (e.g., Mayor||l976) [Nordstrom et ar]|2004t 


[Balser et al.|201 l[|Bovy et al.|2014| Hayden et al.|2014 1 . Metal 
licities of Zsmc = 0.14Zo and Zlmc = 0.43 Zq are assumed for 
the SMC and LMC stars, respectively (Meynet & Maeder 2005 [ 


[Eldridge & Vink|2006||Brott et al.|201 i[ Georgy et al.|2013|l. 


In Fig. 10 we also show a linear regression to the data for 


the full set of analyzed stars, resulting in a power law with an 
exponent of 1.2 ± 0.1. This exponent, however, strongly depends 
on the assumed metallicities. If we alternatively assume, e.g., 
a metallicity of 0.57 Zq for the LMC and 0.29 Zq for the SMC 
( [Meynet & Maeder|2005| l but keep the Galaxy at Zq = 0.014, we 
would obtain a considerably steeper relation with an exponent of 
1.8. An exponent of only 0.9 is achieved for a Galactic and M 31 
metallicity of 0.02, while keeping the LMC and SMC metallicity 
as in Fig.p^ 

For comparison, [Nugis et al.[ ([2007|l suggest an exponent 


on the order of unity, while the studies performed by Crowther 
( [2006| l resulted in an exponent of 0.8 + 0.2. Our favored expo¬ 
nent of a = 1.2 is somewhat higher than the latter value. jWnk &[ 
de Koter[ ( 2005|l theoretically predicted an exponent of 0.86 for 


cool WN stars in the metallicity range 10 ^ < Z/Zq < 1. 

Hydrodynamical stellar atmosphere models (jGrafener & 


Hamann|2005[ [2008[ ) and Monte Carlo wind models fVink & de 


Koter 


2005[l argue that the wind of WR stars is mainly driven by 


radiation pressure through the vast amount of line transitions of 
the iron group elements. This is valid as long as the metallicity is 
not too low. For extremely low metallicities, the CNO elements 
might become more important for the wind driving. In Fig. [m 
we plot the average mass-loss rates in the three galaxies over the 
iron abundances given in Table[T] which is considered to be a 
proxy for the whole iron group. A linear regression to the data 
for the full set of analyzed stars in each galaxy results in a power 
law with an exponent of 1.4 + 0.2. This relation is steeper than 
the M-Z-relation obtained above because the gradient of the iron 
abundance is shallower than the metallicity gradient. Since the 
iron group elements are directly connected to the wind driving 
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Table 4. Coefficients for the modified wind momentum luminosity re¬ 
lations 



log Do 

m 

MW 

0.12 + 0.04 

0.77 + 0.11 

LMC 

-0.34 + 0.04 

1.03 + 0.11 

SMC 

-0.98 + 0.10 

0.26 + 0.59 


mechanism of massive stars, the M-Fe-relation is in principle 
the more appreciable dependency. 

In addition to the direct method used above, we also derived 
the dependence of the mass-loss rate on the metallicity by com¬ 
paring the individual WLRs for the WN stars in the MW, LMC, 
and SMC. In principle, this approach is advantageous because it 
implicitly corrects for the dependence of the mass-loss rate on 
the luminosity. Unfortunately, the WLRs for the WN stars in the 
MW, LMC, and SMC cannot be determined with high accuracy, 
primarily because of the large scatter in the derived mass-loss 
rates ( jHamann et al.[|2006[ [Hainich et al.[|2OT4l l. As a conse¬ 
quence, the slopes of the different WLRs deviate quite strongly. 
Moreover, the exponent of the SMC WLR is very uncertain ow¬ 
ing to the low number statistic. To the samples of each galaxy, 
we fitted a relation of the form 


log£»n 


log Do -H 29.7 + mlog(L/Lo - 5.8) . 


( 6 ) 


The corresponding coefficients can be found in Table]^ while the 
relations, together with the positions of the individual objects, 
are plotted in Fig.[T2] 

Since the modified wind momentum is also a function of the 
terminal wind velocity, the metallicity dependence of Voa needs 
to be considered in the evaluation of the mass-loss rate metallic¬ 
ity relation. For WR stars, however, this dependence is weak, if 
present at all. For example, Niedzielski et al. (20041 found only a 


marginal difference in the terminal velocity between their Galac¬ 
tic and LMC WN sample, studying UV resonance lines such as 
the Heii/l 1640 line. Therefore, we do not consider such a de¬ 
pendence here. Thus, combining Eqs. (j^ and (j^ we obtain 




ccZ" 


(7) 
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Fig. 12. Modified wind momentum versus luminosity for the WN pig. 13. The A log values plotted over the difference of the loga- 
stars in the Galaxy, LMC, and SMC. The dashed lines represent linear nthm of the metallicity. The dashed line is a linear regression, 
regressions to the respective samples. 


while, the slope of the mass-loss rate metallicity relation can be 
written as 


logD 


logZ 


mom logDo + mlog(L/Lo) 

-I- const. = - :-::-h const. . (8) 


logZ 


As the exponents of the WLRs are not identical, one 
has to evaluate Eq. ([^ for a specific luminosity. We choose 
log(L/Lo) = 5.8, which is right in the middle of the overlap¬ 
ping range of the different WLRs. In Fig.[T^ we plot the differ¬ 
ences between the log Dmom values for the three WLRs versus 
the corresponding differences between the logZ values. Since 
AlogDtnom vanishes for AlogZ = 0, the slope of the resulting 
line fit gives the exponent of the mass-loss rate metallicity rela¬ 
tion. With a - 1.3 ± 0.2, it is close to the value directly derived 
from the relation shown in Fig. 10 We note that the quoted er- 


Le = 


KpL 


(9) 


47rcGM 

where is the Thomson opacity. Assuming a completely ion¬ 
ized plasma, which is a reasonable assumption in the deep layers 
of the stellar atmosphere, Le can be written in the form 

L/Lo 


Le^ 10-4-^'3(1 +Zh) 


M/Mo 


( 10 ) 


ror only characterizes the quality of the fit, while several addi¬ 
tional uncertainties would come on top. The largest contribution 
follows from the uncertainty in the derived slopes of the WLRs, 
which in turn relate to the adopted luminosity at which the WLRs 
are evaluated. At a luminosity of log(L/Lo) = 5.6 the exponent 
would be a = 1.1 + 0.3, while it is a = 1.5 + 0.5 for a luminosity 
of log(L/Lo) = 6.0. 

5.3.3. Dependence on the Eddington factor 

The strong mass loss and the spectral appearance of WR stars are 
associated with their proximity to the Eddington limit ( |Grafener] 
& Hamann|2008t Vink et al.|2011| l. This proximity can be char¬ 
acterized with the Eddington factor 


Since Xu and the luminosity L are known from the normalized 
line fit and SED fit, respectively, the only unknown quantity is 
the stellar mass. Unfortunately, a reliable estimate of the surface 
gravity and consequently of the stellar mass from the absorption 
lines in the spectra of the SMC WN stars was unsuccessful be¬ 
cause of the low spectral resolution of most optical spectra at 
hand and the strong impact of the stellar wind on these lines. 
Following the procedure applied by Grafener et al.| (2011 1 and 
Bestenlehner et al. ( 2014[ ), a rough estimate of the stellar mass 
of our sample stars (as a function of their luminosity and surface 
hydrogen abundance) may be obtained by means of the mass- 
luminosity relation derived by Grafener et al. ( |201 1[ Eq. 11). In 
Sect.|^ we compare these masses, which in a strict sense are only 
valid for chemically homogeneous stars with solar metallicities, 
with masses obtained from stellar evolution models. 

It is customary to correlate mass-loss rates of WR stars with 
their Eddington factors Ee, which account for electron scattering 
alone. However, scattering of electrons contributes only a frac¬ 
tion to the total radiative pressure exerted on the atmosphere; in 
hot stars with dense winds, the contributions of line and contin¬ 
uum transitions are far from being negligible. With PoWR, we 
are able to accurately calculate the full radiation pressure and the 
corresponding Eddington factor, which is denoted with Erad and 
which generally varies with depth more extremely than Ee. In 
Fig. 14 we plot the mass-loss rates over the Eddington factor Eg, 
and alternatively over the complete Eddington factor a veraged 


2015 I. While 


over the hydrostatic domain, Ejad ( jSander et al. 

Eg values are well below unity, the E^d values are considerably 
higher, pushing these stars significantly closer to the Eddington 
limit already in the quasi-hydrostatic layers of the atmosphere. 
We note that because of incomplete opacities in the model cal¬ 
culations the true Trad values are probably still underestimated. 

The largest uncertainty in the evaluation of Eg and the Erad, 
however, is caused by the uncertain stellar mass. Evolution 
tracks predict that WR stars with surface hydrogen abundances 
below approximately 40 % (by mass) are already core helium- 
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Fig. 14. Mass loss rates of the WN stars in the SMC versus Edding¬ 
ton factor. The abscissa has two alternative meanings: the standard Ed¬ 
dington factor Ee which only accounts for the radiation pressure on free 
electrons (stars represented by diamonds), and the Eddington factor Tjad 
which includes the total radiation pressure on all absorbers, averaged 
over the quasi-hydrostatic layers (asterisks). 


Fig. 15. The mass-loss rate of the WN stars in the SMC vs. the Edding¬ 
ton factor. The meaning of the colored symbols is the sa me as in Fig[7] 
The open symbols refer to the OB-type stars studied by |Mokiem et al.| 
P007[ l. The gray filled symbols depict the position of the SMC WN 
stars (with Xh < 0.5) when the mass is calculated based on a mass- 
luminosity relation for hydrogen free stars. The dashed line is a linear 
fit to the OB-type stars and the WN stars (colored symbols). 


burning ( |Meynet & Maeder|2005[ [Georgy et ^ 2013| ). However, 
these tracks do not reproduce the HRD positions of our stars 
(see Sect.j^ and can therefore not be used to extract a mass- 
luminosity relation. As an extreme assumption, one may adopt 
the M-L-rel ation for hydrogen-free stars provided by |Grafener] 


et al. 


(20111. However, this leads to Trad values that exceed unity 


and do not allow for the existence of quasi-hydrostatic layers. 
Alternatively, our program stars might be explained by homoge¬ 
nous evolution (see Sect.j^. We therefore choose to calculate 
the stellar masses using the mass-luminosity relation for homo¬ 
geneous stars provided by|Grafener et al. (|201 l|l. 


In Fig. 15 we plot the mass-loss rates versus the Eddington 
factors for our sample stars as well as for the OB-type stars 
studied by Mokiem et al. ( 2007) 1. For consistency, the masses 
of the OB-type stars were estimated with the same method as 
those of the WN stars. In comparison to the OB-type sample 
presented by Mokiem et al. (2007 1 , the SMC WN stars exhibit 
about the same mass-loss rates at the equivalent Fe values. The 
fit to the combined sample of OB-type and WN stars results in 
a power law with an exponent of 1.9 + 0.4. In this context, the 
WN stars in the SMC seem to constitute an extension to the M- 
Fe-relation observed for the OB-type stars. Thus, the strong de¬ 
pendence of the WN mass-loss rate on Fe predicted by Grafener 


|& Hamamil ( |2008| and |Vink et al.| ( [2011| l is not observed in the 
SMC. However, because of low number statistics, these state¬ 
ments are somewhat speculative. 

According to 


Vink et al. 


(20111 a steepening of the M-Fe 


relation is expected for WR stars in comparison to OB-type stars 
because the former objects are in the multi-scattering domain. 
The first observational indications of a steepening of the M-Fg- 
relation were found by |Bestenlehner et ah] ( |2014 1 for stars in 
30Doradus. However, this trend is not seen in the SMC, which is 
consistent with the fact all SMC WN stars fall below the single¬ 
scattering limit. 


For those objects in our sample with hydrogen abundances 
below Ah = 0.5 (mass fraction), we also include in Fig. 15 the Fg 
values obtained with the mass-luminosity relation for hydrogen- 
free stars, indicating the range of possible Eddington factors. For 
the latter Fg values, the observed mass-loss rates of our program 
stars are even lower than the tentative M-Fg-relation (dashed line 
in Fig.[T5]l. 

5 . 4 . WN mass-loss prescription 

The large number of WN stars analyzed in the MW, M 31, the 
EMC, and the SMC should provide the bases for a statistical 
fit to the mass-loss rate as function of fundamental parameters. 
We adopt a mass-loss prescription that is a function of the stel¬ 
lar luminosity L, the stellar temperature T,, the surface helium 
abundance X^e, and the metallicity Z in the form of a power-law 


log(M) = Cl 


Cl log(L) -H C3log(7’.) 
C 4 log(AHg) + C 5 log(Z). 


( 11 ) 


Based on datasets derived from spectral analyses, the free pa¬ 
rameters (Cl - C 5 ) are constrained by means of The best 

fitting parameters are summarized in Table]^ When a specific 
parameter is not considered in the fit, it is depicted by a dash. 

In Table]^ the first line gives the coefficients obtained for the 
whole set of analyzed WN stars, comprising 153 putatively sin¬ 
gle WN stars from Hamann et al. ( 2006|l, Martins et al. ( 2008)1, 
jLiermann et al. ()2010)), jOskinova et al.)^013) ), jSander et ah) 
\20\A) , and Hainich et al. ( 2014 1. The mean deviation between 
the fit and the individual measurements is 0.27 dex and thus rel¬ 
atively large. 

If we restrict the fit to the hydrogen-free objects (Ah < 0.05, 
second line in Table]^ the scatter is comparable (cr = 0.26 dex). 
For the sub-sample of the hydrogen-rich objects (third line), 
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Fig. 16. The mass-loss rate of the WN stars in the SMC, LMC, M31, 
and MW vs. the luminosity. The meaning of the symbols is the same 
as in FiglT] Highlighted are the hydrogen free as well as the stars with 
large surfece hydrogen abundances. 


Table 5. Coefficients obtained by a ;t'^-fit of Eq. 1 11 i to the stellar pa¬ 
rameters of the complete sample of 153 WN stars, and of the sub¬ 
samples of the hydrogen-free (56 stars) and the hydrogen-rich objects 
(14 stars), respectively. 


Ah 

Cl 

Cj 

Cl 

C4 

Cs 

0-0.7 

-5.13 

0.63 

-0.23 

1.3 

1.02 

<0.05 

-7.15 

0.90 

-0.13 

- 

0.97 

>0.5 

-6.17 

0.80 

-0.44 

- 

0.75 


however, the scatter is significantly smaller (cr = 0.07 dex). Both 
these fits exhibit similar M-L relations, with exponents in the 
range of 0.8-0.9. In comparison to the mass-loss relation derived 
by|Nugis & Lamers] ( |2000| l, |Hamann et al.| ( [l995| l, and[H 


et al. (|2014[), all of the prescriptions given in Table|^ posses 




shallower dependence on the luminosity. Moreover, the relation 
for all objects also exhibits a weaker dependence on the helium 


abundance (M oc compared to M oc obtained 


2.22 , 


lby|Nugis 1 


|& Lamers]|2000| l. While the temperature dependence is almost 
negligible in the hydrogen-free case, it is slightly negative for 
the hydrogen-rich stars. Interestingly, the exponent obtained for 
the M-Z relation is smaller for the latter objects. In each case, 
however, the metallicity dependence obtained from the multidi¬ 
mensional fits is shallower than the one-parameter fit conducted 
in Sect. l5.3.2l 

Theoretical mass-loss prescriptions for WN stars were de¬ 
rived by [Grafener & Hamann| ( |2008| l and |Vink et ah] ( |201 l| l. 
Both studies encountered a strong dependence of the mass-loss 
rate on the Eddigton factor. Similar results have been obtained 
by [Grafener et al.| ( |2011| l and |Bestenlehner et al.| ( |20i4) l for WN 
stars in the Arches cluster and 30 Doradus, respectively. There¬ 
fore, we also explored whether this behavior can be observed 
considering a large number of diverse WN stars. As described 
in Sect. 5.3.3| we estimated the mass of the WN stars by means 
of a mass-luminosity relation predicted by|Grafener et al.|(|201 1[ 


their equations 11 and 13) for chemically homogeneous stars and 
solar metallicity. With theses masses, we computed the classi¬ 
cal Eddington factor by means of Eq. ( [T0| ). Incorporating these 
values in the determination of the mass-loss prescription results 
in a relatively weak, or even vanishing, dependence on the Ed¬ 
digton Gamma when also an explicit hydrogen dependence is 
considered in the fit. In each case, the exponent obtained for the 
power-law dependence of the mass-loss rate on Eg is consider¬ 


ably weaker than the values predicted by Grafener & Hamann 
( |2^ and IVink et aL] ( |20TT] l. 


5.5. Stellar rotation 

The rotational velocity of OB stars can be directly determined 
from the profiles of their photospheric absorption lines by means 
of Eourier analyses or spectral fits. Eor the latter, synthetic model 
spectra are convolved with a rotational profile. This method, 
however, is not applicable to the emission line spectra of WR 
stars, since the spectral lines of these objects are formed in the 
stellar wind. Line formation hence requires an integration of the 
radiative transfer equation in three dimensions. A corresponding 


extension to the PoWR code has been introduced by Shenar et al. 

( |20T^ . 

Assuming conservation of angular momentum, the azimuthal 
component of the velocity drops rapidly with increasing distance 
from the star. Adopting such a kinematic model, rotation might 


have significant effects on emission line profiles (Hillier & Miller 
1998[ [Shenar et al.|2014 1. More drastic rotational broadening is 
obtained when assuming co-rotation of the wind up to some dis¬ 
tance. Such co-rotation could be enforced by sufficiently strong 
magnetic fields. 

Since our program stars do not show evidence for wind co¬ 
rotation, we adopt angular momentum conservation and perform 
test calculations with different rotational velocities ysini (de¬ 
fined at R»). Eigure[T7] depicts the obtained line profiles for two 
representative N v lines of SMC AB 12. Erom such comparisons 
with the observed profiles, rotation velocities higher than 200 
km/s can be excluded for all stars of our sample. Sharper lim¬ 
its can be drawn for ABl and AB4 (100 km/s) and for AB2 
(50km/s). The sharper limits rely on characteristic Of lines such 
as N IV T 4060 and N ni A 4640, present in late-type WN spectra. 
The Of lines are usually formed very close to the photosphere 
of the star and their line profiles are thus more strongly affected 
by rotation than lines formed far out in the wind. Therefore, the 
constrains on the rotational velocity are more restrictive for late 
spectral types. [Martins et al.[ ( [2009 j l excluded rotational veloci¬ 
ties much larger than 50km/s for the SMC WN stars by means 
of synthetic emission line spectra convolved with rotational pro¬ 
files. In contrast, our elaborated integration scheme results in a 
high upper limit of about 200 km/s for the rotational velocity of 
the WN 3 stars. 


6. Discussion of the evoiutionary status 

The evolution of massive stars, especially in their late phases, is 
still not completely understood (e.g.,[Langer et al.|1994[[Vanbev 


eren et al.|1998t Langer|2012 Sander et al.|2014j Meynet et al. 

2015[ l. Comparing empirical results with theoretical predictions 
is the basic instrument to unveil the evolutionary status connec¬ 
tions. While stellar evolution models achieved progress in repro¬ 
ducing the observed ratio between the numbers of WN to WC 


stars ([Meynet & Maeder[[2003[ [Eldridge & Vink[[2006] Georgy 

et al.| 2012 | l, they still partly fail to reproduce the range of lumi¬ 

nosities and temperatures observed for the WR population in the 


Article number, page 13 of 30 










































































A&A proofs: manuscript no. paper 



4580 4600 4620 

X,[A] 



4920 4940 4960 

?.[A] 


Fig. 17. Comparison of synthetic spectra for different rotational veloci¬ 
ties, accounting for the full 3D line formation in the wind under the as¬ 
sumption of angular momentum conservation, v sin i refers to the photo¬ 
sphere. The left panel shows the N v d/l 4604,4620 lines of SMC AB 12, 
while the right panel depicts the N v T/i 4934,4944,4951 lines of the 
same object. The blue full line is the observation. The red dashed line, 
green dotted line, and black dashed-dotted line correspond to a rota¬ 
tional velocity of nsini = Okm/s, 200km/s, and 300km/s, respec¬ 
tively. 
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Fig. 18. The HRD of the single WN stars compared to the recent stellar 
evolution tracks (with rotation) calculated by the Geneva group ( |Georgy| 
|et al.|2013^ for a metallicity of 0.14 Zg. The color coding of the symbols 
and the tracks during the WN phase reflects the hydrogen abundance at 
the stellar surface (see inlet). The labels close to the tracks refer to the 
initial mass. 


LMC ( Hainich et al.|2014) l and the Galaxy ( Hamann et al.|2006 
[Sander et al.|2012[ ). Moreover, the current stellar evolution mod¬ 
els predict a considerable number of WR stars in a parameter 
range where no object is actually observed (Hainich et al.|2014 


Hamann et al.|2006| l. In the following, we test whether the same 
holds for the putatively single WN stars in the SMC. 

Figure[T^ shows the HRD position of the WN stars in the 
SMC in comparison to single star evolution tracks published re¬ 
cently by the Geneva group ( Georgy et al.||20l'^ . These tracks 
are calculated for an initial metallicity of Z = O.ITZq. Based 
on the stellar temperature and the surface chemical composi¬ 
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Fig. 19. Same as Fig.jlSl b ut with evolution tracks for Z = 0.28 Zq from 
Meynet & Maeder (2005 ' which also account for rotation|^ 


tion, the discrete WR phases are distinguished by different col¬ 
ors and line styles. The pre-WR phase is defined by a hydrogen 
mass-fraction in excess of Xh = 0.55 (thin black lines). At hy¬ 
drogen abundances below Xu - 0.55 and stellar temperatures 
above T* = 25 kK, the track is considered to represent the WNL 
stage (bold red line). Once the hydrogen abundances drops be¬ 
low Xu - 0.05, we classify the star as a WNE (thick green line). 

The lowest initial mass for which a track reaches the WR 
phase is 85 Mq. Interestingly, the tracks for 85 Mq and 120 Mq 
enter the WN stage for the first time already at the end of the 
main sequence. Subsequently, both tracks evolve toward cooler 
surface temperatures in the regime of the supergiants before they 
experience another blueward evolution and enter the WN stage 
for a second time. While the second WN phase is only short 
for the 85 Mq model, the 120Mo model appears as WNE even 
during carbon burning, close to the end of the evolution track. 

In contrast to these predictions, the SMC WN stars are found 
in a region of the HRD (log(L/Lo) = 5.6...6.1) that corre¬ 
sponds to considerably lower initial masses in the range from 
about 30 Mq to 60 Mq. The comparison illustrated in Eig. 18 
clearly elucidates that the parameter regime occupied by the 
SMC WN stars is not reproduced by the theoretical prediction 
of the Geneva evolution models for single stars with an initial 
metallicity of Z = O.ITZq. 

In Eig.fT^ we compare our sample with the older generation 


of the Geneva evolution models computed by Meynet & Maeder 
( |2005| l, which also account for the effect of rotation but are cal- 

^ The differences between the tracks shown in this figure and the tracks 
plotted in the left panel of Fig. 6 in jMartins et al.| ( |2009^ are attributable 
to the different kinds of effective temperatures given by [Meynet ^ 
[Maeder[ l [2005^ . Martins et al.[ ( [200^ used the temperature that is cor¬ 
rected for the thickness of the WR wind, whereas we plot the tempera¬ 
ture that refers to the hydrostatic radius. The latter temperature is com¬ 
parable to the stellar temperature T, used in stellar atmosphere codes 
like PoWR or CMFGEN. 
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Fig. 20. Same as Fig. |l8[ but with non-r otating evolution tracks for Z = 
0.28 Zq calculated by |Eldridge & Vink| ( |2006| l. 


culated for a higher initial metallicitiy of 0.28 Zq. In contrast to 
the recent Geneva models for Z = 0. 14 Zq ( |Georgy et al.|2013| l, 
the tracks plotted in Fig.[^ evolve considerably further to the 
blue part of the HRD. Some of these tracks even reach the WC 
phase (gray line), defined when the carbon mass-fraction at the 
surface exceeds Zc = 0.2. Moreover, compared to the evolution 
models published by Georgy et al.| ( |2013| ), the tracks drawn in 
Fig.[^ (|Meynet & Maed^ 2005[) predict a significantly lower 


initial mass limit for reaching the WR phase of 40 Mq. However, 
the agreement between our empirical HRD positions and the 
models from |Meynet & Maeder] ( |2005| l is still relatively small. 
Only the two most luminous stars in our sample (SMC AB 9 and 
1) are reproduced. The relatively low luminosities of SMC AB 2 
and 10 suggested that these stars evolved from an initial mass 
of about 30 Mq, while the corresponding track never turns into a 
WR star. 

In the above comparison remains unclear whether the 
adopted metallicities are causing the different results, or whether 
other details of the model calculations are the cause. We also 
compared the HRD positions of our program stars with the non¬ 
rotating stellar evolution models from Eldridge & Vink| (2006 


see Fig.[20)l, which are calculated for an initial metallicity of 
0.28 Zq. The mass-loss rates which they adopted for their non¬ 
rotating models in the pre-WR phases is comparable to the one 
used by |Meynet & Maeder| ( |2005] ) for their rotating models. Dur¬ 
ing the WR phase, the models described in [Eldridge & Vink| 
( 12006| l account for a metallicity scaling of the mass-loss rate, 
while such a dependency was not considered in the older gen¬ 
eration of the Geneva models ( [Meynet & Maed^|2005[ ). Tak¬ 
ing into account that the total WR lifetime is nearly independent 
from the metallicity scaling of the mass-loss rate during the WR 
phase ( [Eldridge & Vink [2006[l, the comparison of these tracks 
with the ones shown in Figs. [18|and[T^ might give some indica¬ 
tion of the impact of the metallicity. 


For the stars in our sample with a luminosity of log (L/Lq) > 
5.8, there is a good agreement with respect to the tracks 
from Eldridge & Vink ( 2006jl, whereas the two faintest objects 


(SMC AB2 and 10) cannot be reproduced at all. Although the 
tracks from [Eldridge & Vink[ ( [2006| ), compared to the Geneva 
tracks, coincide better with our results, these models also require 
a minimum initial mass of 45 Mq to reach the WR phase and 
therefore do not reach the WN phase for luminosities as low as 
those of SMC AB 2 and 10. 

The higher the metallicity of the stellar evolution models, 
the higher the degree of agreement that can be achieved with 
the HRD positions of the WN stars in the SMC. Our sample is 
best reproduced with tracks for Z = 0.56 Zq from Meynet & 
[Maeder[ ( [2005 [ l, which account for the implications of rotation. 
The model atmospheres used in our analysis have a metallicity 
of Z = 0.13 Zq. Our relatively young WN stars might have been 
born with a slightly higher metallicity than this average value 
for the SMC. However, we do not see evidence for a consid¬ 
erably higher iron abundance than the value given in Table[2 
An exception is SMCAB4, which might be reproduced with 
stellar atmosphere models of higher metallicities as well. The 
enhanced nitrogen abundance observed (see Table [A.3l ) can be 
well understood in terms of self-enrichment with products from 
central helium burning. In conclusion, we suggest that it is not 
the metallicity itself, but some other effect which is not prop¬ 
erly accounted for in the evolution models that are causing the 
difference between the tracks and the empirical HRD positions. 

Apart from the discrepancies described above, another dis¬ 
agreement between our analysis and the stellar evolution models 
is the surface hydrogen abundances. When the evolution tracks 
cross the empirical HRD positions of our SMC WN stars, they 
predict a considerably lower hydrogen abundance at the surface 
than observed. This coincides with the fact that the mass-loss 
rates derived in our analysis (see Sect. HU are on average lower 
than the values given by the mass-loss relation from [Nugis & 


[Lamers| ( |2000[ l, which is often applied in stellar evolution mod- 
elsj^Thus, the discrepancy in the hydrogen abundance might be 
caused by too high WR mass-loss rates in the evolution models, 
which remove too much of the hydrogen-rich envelope. There¬ 
fore, higher mass-loss rates in the WR phase would not bring 
the stellar evolution tracks into accordance with the empirical 
HRD positions of the SMC WN stars. We conclude that either 
the mass-loss in a pre-WR stage is significantly higher (via either 
continous mass-loss or episodic bursts), as suggested by e.g.. 


ISander et al. 

(|2012|l to explain the disagreement between WC 

stars and evo 

ution models at Galactic metallicity (see also[Van-[ 

beveren et al. 

1998 2007[ Georgy 2012|l, or that standard models 


for single star evolution still miss essential physical ingredients. 

Such a missing ingredient could be very fast rotation, which 
lead to full internal mixing and quasi-homogeneous evolution. 
This type of evolution proceeds roughly speaking directly from 
the hydrogen zero age main sequence to the helium main se¬ 
quence (e.g., Maeder|1987[[Langer|1992[[Yoon & Langer|2005[ 


Woosley & Heger|2006| l. Homogeneous evolution is expected to 
be more common at low metallicities, since rotational mixing is 
more efficient ( [Brott et al.|2011[ ) and angular momentum loss is 
lower because of the reduced wind mass-loss (see Sect. [5.3.2] i. 
Thus, rapid initial rotation will be conserved for a longer period, 
facilitating the evolution toward the blue part of the HRD even 
at relatively low initial masses. 


In some cases the more recent description derived b y|Grafener & 
Hamannl (|2008b is also used (e.g., [Ekstrdm et al.|2012| [Georgy et al. 
2013) . 
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Fig. 21. The HRD of the single WN stars compared to the evolution 
tracks that account for very rapid initial rotation > 500km/s, Brott| 
|et aL|20TT| l. These evolution models assume an initial metallicity of 
0.14Zo and were calculated roughly until the central hydrogen exhaus¬ 
tion. The labels close to the tracks refer to the initial mass. The color 
coding in the same as in Fig.[T^and Fig.[T^(see also the inlet). 


For SMC AB 1 and 2, this scenario was already proposed by 
[Martins et al.| ( |2009) , mainly because of the high hydrogen abun¬ 
dance and the location of these two objects in the HRD. As ev¬ 
ident from Fig.|^ stellar evolution tracks ( |Brott et'aL |201 l| l, 
when accounting for vary rapid initial rotation (njot > 500km/s) 
and a metallicity of 0.14Zo, can at first glance explain the HRD 
positions of most SMC WN stars. Taking a closer look at the 
evolution of chemically homogeneous stars, however, reveals 
that the hydrogen abundance of these models in the correspond¬ 
ing parameter range is considerably lower than derived from 
the spectral analysis. The 50 Mq track close to the position of 
SMC AB 1, e.g., predicts a hydrogen mass-fraction below 9%, 
while our analysis (in agreement with the hndings of |Martin^ 


et al. 20091 results in a value of Xu = 50 %. For comparison, 
the 80 Mq track from Eldridge & Vink ( 2006[ ) exhibit a hydro¬ 
gen abundance of about 15 % close to the position of SMC AB 1. 
The disagree ment between our results and the evolution tracks 
published by |Brott et al. (20111 is conhrmed by the BONNSA^ ^ 
Bayesian tool ( [Schneider et ak|2014[ ), which allows us to qua' 
itatively match the Bonn stellar evolution models (jBrott et'a 


[201 1[ [Kohler et al.[[2015| ) to the results obtained from our ste 
lar atmosphere analyses, accurately accounting for observational 
uncertainties. With the exception of SMC AB4, the parameters 
derived in our analysis could not be consistently reproduced by 
this tool. Thus, the discrepancy between the hydrogen abundance 
predicted by the evolution models and the values derived by 
means of the spectral analysis does not support the scenario of 
homogeneous evolution. 


^ The BONNSAI web-service is available at www.astro.uni- 
bonn.de/stars/bonnsai. 
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Fig. 22. The evolution of the equatorial rotational velocity as a function 
of time for the stellar evolution models calculated by [Brott et al.[j2011[ 
solid lines) and a 60 Mq track (dashed line) published by [Meynet &[ 
[Maeder[ ([2005 |l. Both studies account for very rapid rotation (Urot.ini > 
500km/s) and an initial metallicity of 0.14 Zq. The labels close to the 
tracks refer to the initial mass. The color coding is the same as in 
Figs.[^and[T^(see also the inlet). 


From homogeneous models, it is expected that the rotational 
velocity is substantial (urot > lOOkm/s) during a considerable 
part of the stellar evolution. For the quasi-homogeneous mod¬ 
els shown in Fig.|^ ( Brott et ar][201 Ij yrot.ini > 500km/s), the 
evolution of the equatorial velocity as a function of time is illus¬ 
trated in Fig.|^ The rotational velocity in these models never 
falls below lOOkm/s. For comparison we include in this hgure 
a 60 Mq track published by Meynet & Maeder (20051 that also 
accounts for fast rotation (uiot.ini = 500 km/s) and quasi homoge¬ 
neous evolution. In contrast to the tracks calculated by [Brott etal.[ 
( 2011| l, the Geneva track exhibits considerably lower rotational 
velocities during the WR phase. As discussed in Sect. [5.5[ rota¬ 
tional velocities as high as 200 km/s at the stellar photosphere 
could not be excluded for the analyzed WN 3 stars, as long as 
the winds do not co-rotate with the photosphere. Thus, homoge¬ 
neous evolution remains a promising possibility to explain the 
origin of the WN stars in the SMC. 

The influence of binary evolution on the rotational velocities 
of massive stars was recently studied by |de Mink et al. ( 2013| l by 
means of a rapid binary evolution code. Assuming a uniform ini- 
tial distribution of the rotational velocity with iirot < 200 km/s, 
these authors show that the rapid rotators observed in popula¬ 
tions of main-sequence stars might consist exclusively of stars 
that already underwent Roche-lobe overflow. Hence, initially 
fast rotation could be very rare among massive stars. However, 
[Cantiello et al.[ ( [2007| l demonstrate that binary interaction can 
provide high rotation rates and quasi-homogeneous evolution 
very similar to the single star evolution models. 

The effects of binary evolution on WR populations were 
studied by various authors ([Paczyhski[[1967[ [Vanbeveren et al.[ 


2007[ Eldridge et al.[[2008| 2013j l. The initial mass required to 


enter the WR phase is considerably reduced, if angular momen¬ 
tum and mass transfer by Roche lobe overflow occurs in the 
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course of the evolution. Therefore, binary evolution can in prin¬ 
ciple account for the HRD positions observed for the SMC WN 
stars. The current binary status of the WR stars in the SMC was 
carefully investigated by |Foellmi et al.| ( |2003| see also Sect.|^, 
inferring that the stars in our sample are not short periodic bi¬ 
naries (with a lower certainty for SMC AB 9). However, accord¬ 
ing to |de Mink et al.| ( |201 1[ |2014| l, a significant fraction of the 
stars that are now single might be products of previous binary 
interaction. Thus, our sample of apparently single stars could be 
polluted or even completely consist of objects affected by binary 
evolution. A study of the properties of the SMC WR binaries is 
currently underway, and should shed light on the systematic dif¬ 
ferences between the single and binary WR population (Shenar 
et al., in prep.). 


7. Summary and Conclusions 

We performed a comprehensive study on the putatively single 
WN stars in the SMC. Optical and, if available, UV spectra were 
analyzed by means of PoWR model atmospheres. The sample 
stars prove to be significantly different from their counterparts 
in the LMC and the Galaxy. In contrast to these galaxies, all 
SMC WN stars show a significant amount of hydrogen at their 
surface (Xh > 0.2). With the exception of SMC AB 4, all objects 
exhibit weaker winds. The mass-loss rates are on average lower 
by about a factor of four and twelve compared to the WN stars 
in the LMC and MW, respectively. This can be interpreted as a 
correlation between the metallicity Z and the mass-loss rates in 
the form M oc Z'-^ This finding substantiates that the winds of 
WN stars are driven by radiation pressure that is intercepted by 
metal lines. We also derived multidimensional fits to express the 
mass-loss rate as a function of the stellar parameters. 

Standard stellar evolution models only partly prove to re¬ 
produce the empirical HRD positions of the WN stars in the 
SMC and their observed surface hydrogen abundance. In gen¬ 
eral, these evolution tracks require an initial mass to reach the 
WR phase that is relatively high, so that they cannot reproduce 
the low luminosities obtained for some of the analyzed stars. Our 
results actually agree better with evolution models calculated 
for higher metallicity. This could indicate enhanced or episodic 
mass-loss in a pre WR stage. Higher mass-loss rates in the WN 
phase do not agree with our results, since the mass-loss rates ap¬ 
plied in stellar evolution models are already higher than we find 
empirically. Moreover, the high hydrogen abundances derived in 
our analysis are in general not reproduced by the stellar evolu¬ 
tion models. Enhanced mass loss in an earlier phase of the stellar 
evolution, e.g., in a luminous blue variable (LBV) phase, might 
explain the blueward evolution of our sample stars, while it could 
preserve enough of the hydrogen-rich atmosphere to explain the 
observed surface abundances. 

An alternative explanation for the evolutionary origin of 
the SMC WN stars is by fast rotation which entails a quasi- 
homogeneous evolution. Corresponding tracks can reproduce 
the observed HDR positions. However, lower mass-loss rates in 
the quasi homogeneous evolution models are necessary to bring 
the predicted hydrogen abundance in accordance with the obser¬ 
vations. 

We also cannot completely exclude that our sample stars 
have been affected by binary evolution. With respect to the high 
hydrogen abundances of the SMC WN stars two scenarios seem 
possible: Either the objects that we currently observe as single 
WN stars were mass gainers in interacting binary systems, which 
were already disrupted, or these objects were formed via merger 
events. 


Summarizing, the Wolf-Rayet population in the Small Mag¬ 
ellanic Cloud is found to be strikingly different from correspond¬ 
ing stars in environments with higher metallicity (our Galaxy, the 
LMC and M31). This reveals fundamental differences in how 
massive stars evolve under metal-poor conditions. With a metal¬ 
licity that is about six times lower than the Sun’s, the Small Mag¬ 
ellanic Cloud comes close to conditions that are typical for the 
early Universe. 
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Appendix A: Additional tables 


Table A.2. The spectra used in the analysis of the putatively single WN 
stars in the SMC 


Table A.l. Atomic model used in the stellar atmosphere calculations 


Ion 

Number of levels 

Number of lines 

Hi 

22 

231 

Hii 

1 

0 

He I 

35 

595 

Hen 

26 

325 

Hem 

1 

0 

Ni 

1 

45 

Nil 

38 

703 

Niii 

36 

630 

Niv 

38 

703 

Nv 

20 

190 

Nvi 

1 

91 

Ci 

1 

105 

Cii 

32 

496 

Cm 

40 

780 

Civ 

21 

300 

Cv 

1 

10 

On 

37 

666 

Om 

33 

528 

Oiv 

25 

300 

Ov 

36 

630 

Ovi 

16 

120 

Ovii 

15 

105 

Piv 

12 

66 

Pv 

11 

55 

Pvi 

1 

0 

Gii“ 

1 

0 

GiiC 

8 

25 

Giv'’ 

11 

49 

Gv“ 

13 

69 

Gvi“ 

17 

121 

G vii^' 

11 

52 

G viiC 

9 

34 

Gix'' 

9 

35 

Gx“ 

1 

0 


Notes. G denotes a generic atom which incorporates the following 
iron group elements Fe, Sc, Ti, Cr, Mn, Co, and Ni. The correspond¬ 
ing ions are treated by means of a superlevel approach (for details see 
IGrafener et al.|200^ . 


SMC AB 

FUSE 

lUE 

UVES 

EMG" 

TM'’ 

1 


X 

X 

X 

X 

2 

X 

X 

X 

X 

X 

4 

X 

X 

X 

X 

X 

9 

X 



X 


10 

X 



X 


11 




X 


12 




X 



Notes. A X i ndicates observation use d in our analysi s. observations 


conducted by Foellmi et al. 

(|2003j and 

Foellmi 

2004 

flux calibrated 

observations conducted by 

Torres-Doc 

gen & V 

assey 

(fT^ 


Table A.3. Abundances derived in the analysis of the putatively single 


WN stars in the SMC (mass fraction in 

percent) 


SMCAB 

An 

2fc 


Ap 

1 

0.225 

- 

- 

- 

2 

0.45 

0.005 

- 

0.00006 

4 

0.675 

- 

- 

0.00006 

9 

0.45 

- 

0.00075 

- 

10 

0.45 

- 

0.0025 

- 

11 

0.225 

- 

- 

- 

12 

0.9 

- 

- 

- 


Notes. When a specific abundance is not determined in the fit, it is de¬ 
picted by a dash. 


Table A.4. Number of ionizing photons and Zanstra temperatures for 
the putatively single WN stars in the SMC 


SMC AB 

log 2 

[s-i] 

Hi 

^Zanstra 

[K] 

He I 

log 2 

[s-i] 

log 2 
[s-i] 

Hen 

Fzanstra 

[K] 

1 

49.94 

87.7 

49.68 

47.68 

70.2 

2 

49.40 

51.5 

48.78 

38.40 

20.3 

4 

49.61 

47.4 

48.94 

37.51 

18.6 

9 

49.90 

106.6 

49.70 

48.49 

96.3 

10 

49.50 

104.9 

49.30 

48.10 

95.8 

11 

49.72 

99.3 

49.49 

47.88 

82.0 

12 

49.73 

118.0 

49.55 

48.57 

109.3 


Table A.5. Mass-loss rates and luminosities derived by iNugis et al.l 
( |2057t 


SMCAB 

log AfjVugis 
[Mo/yr] 

logM 

[Mo/yr] 

log 

[^o] 

logL " 
[Lo] 

1 

-5.41 

-5.58 

5.88 

6.07 

2 

-5.58 

-5.75 

6.24 

5.57 

9 

-5.62 

-5.65 

5.62 

6.05 

10 

-5.7 

-5.64 

5.38 

5.65 

11 

-5.69 

-5.56 

5.59 

5.85 

12 

-5.85 

-5.79 

5.23 

5.90 


Notes. this work 
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Appendix B: Comments on individual stars 


SMC AB1 (alias: AzV 2a, SMC-WRl): We achieved the best fit 


with a stellar temperature of 79.4 kK, while Martins et al. (2009 1 
derived a value of 65.4kK. At this low temperature, however, all 
observed N v lines are considerably underpredicted by our mod¬ 
els. In addition, the synthetic spectra of these models posses a 
distinct N iv 4 4060 emission line that is not observed for this 
object. The N iv emission in the model spectra already vanishes 
at approximately 70 kK. Because of a better fit of the N v lines 
we favor, however, the model with a slightly higher temperature. 
Another indication of the higher temperature is the Civ4 1551 
line, which also vanishes at approximately 70 kK and which is 
indeed not observed for this object. The higher temperature in 
comparison to the findings by [Martins et al.| ( |2009| l also results 
in a stellar luminosity which is 0.32 dex higher. Although both 
studies are based on the same ESO VLT UVES observations, 
the reduced data have some minor differences. Eurthermore, the 
differences in the stellar parameters are partly due to the fact 
that the models do not reproduce all lines in the observed spec¬ 
trum consistently. Hence the derived stellar parameters depend 
on which lines are considered to be most significant for the line 
fit. We note that the mass-loss rate is almost identical in both 
studies. 


SMC AB 2 (alias: AzV 39a, SMC-WR2): |Testor| ( |200T] l found 
this star to be embedded in the HII region N 28. The stellar pa¬ 


rameters of SMC AB 2 were already studied by Martins et al. 
( [2009| l. Within the expected error margins, the stellar temper¬ 
ature and the luminosity obtained in this study are compatible 
with the findings of the previous authors. Martins et al.|P009] l 
assume a clumping factor of D = 10, while our fit points to a 
lower value of about four. Taking this into account, the mass- 
loss rates obtained in both studies deviate by only 0.1 dex, which 
is within the expected error range as well. 


SMCAB4 (alias: AzV 81, Sk41, SMC-WR4): It is the only 
WN6h star in the SMC which appears to be s ingle ([Eoellmi et ah 
2003 . It was analyzed before by Crowther] ( 2000|l and Martins 


et al. (|2009|l. Compared to the findings of fhese two studies. 


our analysis results in a slightly higher stellar temperature (by 
approximately 2-3 kK) and a mass-loss rate that is lower by 
about 0.1 dex when the different clumping factors are taken into 
account. The luminosity deduced from the SED fit lies within 
the range determined by the previous analyses ( |Crowther|2000[ 
[Martins et al.|[2009) l, while the deviation amounts to approxi¬ 
mately 0.1 dex in both cases. We also note that the hydrogen 
abundance deduced in our analysis is lower by 5 % and 10% 
in relation to the results of |Crowther| ( [2000) 1 and [Martins et al.[ 
(2009|l, respectively. 


SMC AB9 (alias: MVD 1, SMC-WR9): This object belongs to 
the hottest stars in the sample. The temperature determination 
is base primarily on the antagonistic behavior of the nitrogen 
lines and the helium lines (especially Hen45412). While the 
nitrogen lines tend to favor lower temperatures, a superior fit 
of the He II 4 5412 line requires higher model temperatures. We 
note that a weak O iii 4 5007 nebular line is visible in the optical 
spectrum, which could indicate an associated nebula. 

Unfortunately, the quality of the fit in the EUSE range is 
poor, since the continuum shape of the EUSE spectra cannot 
be reproduced by the models. This seems to be an issue of the 



Fig. B.l. HST-WFPC2 images of the region within about 8 " distance of 
SMC AB 11. Left panel: Image taken in the F606W filter. Right panel: 
Image taken in the F300W filter. The red circle marks the position of 
SMC AB 11. The images were retrieved from the HST archive. North is 
up and east is left. 


observation rather than a model deficiency. The observed flux 
shortward of 1000 A is approximately a factor of two lower than 
predicted by the model. Since no significant ionization edge is 
expected in this spectral region, we assume that the observation 
at hand is affected by inaccurate flux in at least one of the two 
SiC channels (SiC refers to the silicon carbide coating of the cor¬ 
responding mirror) of the FUSE satellite. This is consistent with 
a bump in the spectrum at about 1085 A, where a gap between 
the segments of the LiF channels (LiF refers to the lithium fluo¬ 
ride coating of the corresponding mirror) is covered by the SiC 
channels. 


SMC AB 10 (alias: SMC-WRIO): The fit of the FUSE spectral 
range suffers from the same caveats as described for SMC AB 9 
(see corresponding comment). In contrast to the other objects in 
our sample, the best fitting model assumes a velocity field in the 
form of a double-y6 law ([Hillier & Miller||1999|l as defined by 
[Todt et al.[ ( p0T5l l with /3i = 1 and /32 = 4. For SMC AB 10, the 
choice of a double-/! law is mainly motivated the exceptional line 
profile of the He ii 4 4686 line that has a broad base and a very 
narrow line peak. With respect to a model using a single-/! law 
with P - I, the mass-loss rate is increased by about 0.1 dex and 
the terminal wind velocity is higher by about 600km/s. 


SMC AB 11 (alias: SMC-WRl 1): A bright red object is located 
about 1.2 " from SMC AB 11 (see also Crowther & Hadfield[ 
2006| l. Since the spectra at hand ( [Eoellmi et al. 2003| l have a 
spatial resolution on the order of 1", we note the possibility of 
a contamination. Figure [BT[ shows two images taken with the 
Wide-Field Planetary Camera 2 (WFPC2) aboard the Hubble 
Space Telescope (HST) in the F300W and F606W filter, respec¬ 
tively. The position of SMC AB 11 is marked with a red circle. 
The HST observations clearly show that the WR star dominates 
the flux in the UV, while both objects exhibit approximately the 
same flux level in the F606W filter. In the IR on the other hand, 
the flux is dominated by the non-WR object. 

The first panel of Figure [CT8| shows the SED fit derived in 
our analysis. Since no flux-calibrated spectra are available for 
this object, the luminosity are derived solely from photometry. 
A variety of photometric measurements is available within a ra¬ 
dius of 2 " from the coordinates given for SMC AB 11 in the 
SIMBAD database. In comparison to the line-of-sight compan¬ 
ion, SMC AB 11 is located farther to the north (see Fig. 13- 
Since the measurements performed by [Zaritsky et al.| ( |2002[ ) and 
[Kato et ar[ ( [2007[ l were obtained at a similar position, we mainly 
used these values as well as the spectrophotometry obtained 
by Crowther & Hadfield ( 2006[ ) to constrain the luminosity of 
SMC AB 11. In Fig.|C.8| the corresponding photometry marks 
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and the model SED for SMC AB 11 are shown by brown boxes 
and a brown dotted line, respectively. Other photometry (Zarit 


sky et al. 2002[ Monet et al.|2003[ Oey et al.|2004[ [Gordon et al. 


2011||Cutri & et al.|2012[ Cutri et al.|2012| l, which presumably 

incorporates contributions from both objects, were used to con¬ 
struct the observed SED that combines both objects (blue boxes 
in Eig. jC.Sj l. In this figure, the red straight line represents the 
model that incorporates both the WN star as well as flux contri¬ 
bution from the non-WR object, which is assumed to be a black 
body (green dashed line). Because of the flux contribution of the 
line-of-sight companion, the uncertainty in the derived luminos¬ 
ity of SMC AB 11 is higher than for the other stars in our sample. 
In the lower panels of Eig. jC.Sj the relative contribution of the 
WN star to the combined spectrum (red dashed line) is shown 
by the brown dotted line. 

Foellmi et al. ( [2003 | l mentions the presence of N iv d 4060 


emission in an ESO-NTT spectrum of SMC AB 11, although nei 
ther the spectrum at hand nor the spectrum published by Massey 
|& Duffy| ( |2001| l exhibits this line. We note that stellar atmosphere 
models with temperatures of about 65 kK would be necessary to 
reproduce this line simultaneously with the N v emission lines. 
This would be considerably lower than the 89 kK derived in this 
analysis. 


SMC AB 12 (alias: SMC-WR12): It is the object with the high¬ 
est surface temperature in our sample. In the same way as the 
other WN3 stars only one ionization stage per element is visible 
in our spectra of SMCAB 12, so that the temperature determi¬ 
nation has to rely on the Nv lines and the Hen lines. At least 
112kK are necessary to reproduce the Heii/15412 line. More¬ 
over at these temperatures the model spectra do not show He ii 
lines originating from the fifth quantum level, which is consis¬ 
tent with the observations. 


Appendix C: Spectral fits 

The following figures show the final spectral fit of each object 
analyzed in this paper. In all plots the upper panel shows the fit of 
the SED, while the lower panels exhibit the fit of the normalized 
optical and UV spectrum, if available. The best-fitting model is 
plotted in red, while all observations are shown in blue. 
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SMC AB 1 WN3ha 

T* = 79 kK logR, = 1.37 
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Fig. C.l. Spectral fit for SMC AB 1 
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SMC AB 2 WN5ha 

T* = 47 kK logR, = 1.63 
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Fig. C.2. Spectral fit for SMC AB 2 
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Fig. C.3. Spectral fit for SMC AB 2 (continued) 
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SMC AB 4 WN6h 

45 kK logRt = 1.30 
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Fig. C.4. Spectral fit for SMC AB 4 
















































































A&A-paper, Online Material p 26 



5500 5750 6000 6250 6500 6750 


X[k] 


Fig. C.5. Spectral fit for SMC AB 4 (continued) 
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SMC AB 9 WN3ha 

T* = 100 kK logRj = 1.15 
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Fig. C.6. Spectral fit for SMC AB 9 
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SMC AB 10 WN3ha 

T* = 100 kK logRj = 1.14 
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Fig. C.7. Spectral fit for SMC AB 10 
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SMC AB 11 WN4h 

T* = 89 kK logR, = 1.35 
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Fig. C.8. The spectral fit for SMC AB 11. The brown boxes in the top panel are the photometry attributed to the WR star alone, while the blue 
boxes refer to photometry that comprises the flux from a companion separated by approximately 1.2" (see Appendix]^ for details). The brown 
dotted line represents the best fitting PoWR model for the WN star, while the green dashed line shows a black body spectrum at 4.6 kK as a rough 
approximation to the nearby red star. The photometry with large apertures (blue boxes) are obviously determined by the red star for A > 8000 A. 
The straight red line is the combined model SED that incorporates the flux from the WN star and the line-of-sight companion. 
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SMCAB12 WN3ha 

T* = 112kK logRj = 1.11 
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Fig. C.9. Spectral fit for SMC AB 12 












